Ph. 



Stochastic evolution of finite level systems: 
classical vs. quantum 



3,4 



D. Chruscinski^* V.I. Man'ko^, G. Marmo^'", and F. Ventriglia 

^Institute of Physics, Faculty of Physics, Astronomy and Informatics, 
■ Nicolaus Copernicus University, Grudziadzka 5, 87-100 Toruh, Poland 

^ I ^P. N. Lebedev Physical Institute, Russian Academy of Sciences, 119991 Moscow, Russia 

D ■ ^Dipartimento di Fisica and MECENAS, 

Universita di Napoli Federico II, 1-80126 Napoli, Italy 



^INFN, Sezione di Napoh, 1-80126 Napoh, Italy 



Abstract 



^ ' Quantum dynamics of the density operator in the framework of a single probabihty vector is 

analyzed. In this framework quantum states define a proper convex quantum subset in an appropriate 
simplex. It is showed that the corresponding dynamical map preserving a quantum subset needs 
not be stochastic contrary to the classical evolution which preserves the entire simplex. Therefore, 
^ ' violation of stochasticity witnesses quantumness of evolution. 

o ' 

m ! 1 Introduction 
^ , 

I The basic idea of quantum tomography is to encode the information of a quantum state in a family of 
. probability distributions. This idea goes back to Pauli [1] who considered the possibility to associate 
] quantum states with marginal probability distributions in the phase space, both for positions and mo- 
menta, as in classical statistical mechanics. In spite that these marginal distributions do not determine 
the quantum state wave function [2] , that idea was realized by considering an infinite family of marginal 
distributions (tomograms). 

The tomographic representation is based on the Radon transform of the Wigner function which relates 
. the measurable optical tomographic probability to reconstruct the Wigner function of the quantum state 
(see [3] for a recent review). Similar tomographic approach may be constructed for finite level quantum 
systems (spin systems) corresponding to finite-dimensional Hilbert space mOIB]. In this case a spin (or 
qudit) density operator p can be described [7] by probability distributions of random spin projection (spin 
tomogram) depending on the direction of the quantization axis. Such family of probability distribution 
encodes the complete information about qudit state and hence a density operator may be uniquely 
recovered out of a spin tomogram. 

Quantum tomography enables one to encode the information about p into the finite family of proba- 
bility distributions (see the next section) and eventually to represent p as a single probability vector living 
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in an appropriate simplex. It is clear that only a convex quantum subset within this simplex correspond 
to legitimate density operators. Using this stochastic representation p — >■ P one may translate a quantum 
evolution of p(t) into a stochastic evolution of P(i). The problem of stochastic evolution and dynamical 
maps of density matrices was studied and reviewed recently in [U [9l [10] . In this way both classical and 
quantum evolution of finite level systems may be described within the same framework. Now comes a 
natural question: suppose that one is given the evolution of a stochastic vector P{t) ~ T(t)P represented 
by a semigroup of maps T{t) satisfying T{t + u) ^ T{t)T{u) for all t,u > 0. Is it possible to discriminate 
wether T(i) describes purely classical evolution of P or a quantum evolution of p encoded in P? It is 
clear that classical dynamics is governed by stochastic matrices T(i). Interestingly, it needs be no longer 
true for quantum dynamics. Some preliminary aspects of this phenomenon were mentioned in 112] 
Note that in the quantum case one requires that T(t) maps probability vectors P within the quantum 
subset of the simplex. In particular it turns out that if the dynamics of p(t) is unitary then T{t) is never 
a stochastic family of maps. This way violation of stochasticity witnesses quantumness of evolution. 

The paper is organized as follows. In Section [2] we review the tomographic approach to quantum 
mechanics. In Section [3] we show how to encode the information about a density operator p into a single 
probability vector P living in an appropriate simplex. Vectors corresponding to legitimate quantum 
states form a convex quantum subset within a simplex. It is shown how the evolution of p(t) induces the 
evolution of P{t). Finally, Section 2] provides several examples of qubit dynamics in the framework of 
6-dinicnsional probability vector P. Final conclusions are collected in the last Section. 

2 Tomography of spin states 

In this section we review the approach [3] to describe the quantum states of a spin (qudit) by means 
of the probability distributions called quantum tomograms. The tomogram W{m\U) corresponding to a 
qudit state with the density operator p is given by the diagonal matrix elements of the unitarily rotated 
density operator in the standard computational basis |m), where —j < m < j and Jzlm) = m|m), with 
Jz being a spin projection operator on z-axis. Thus one has 



where U is a (2j + 1) x {2j + 1) unitary matrix. In the case when the matrix J7 is a matrix of irreducible 
representation of the group SU(2) the tomogram W(TO|n), called a spin tomogram, introduced in |4] [5] 
depends on the variable n = (sin0 cos yj, sinS sin (/s, cos 0). Spherical angles (d,^) determine the point on 
the Poincare sphere. Both tomograms W{m\U) and M^(TO|n) define probability distributions, that is, 
W{m, U) > 0, W{m, n) > 0, and the normalization condition 



These properties mean that the tomograms can be interpreted as conditional probability distributions 
W{m,U) = W{m\U) and W{m,n) = W{m\n) which can be associated with a joint probability dis- 
tributions V{m,U) or V{m,n). The construction of the probability distribution V{m,n) in the form 
V{m,n) = iy(TO, n)n(n), where n(n) is a probability distribution on the sphere 5^, was discussed for 
particular cases in [T3J [H] and the probability vector with components W{m,Uj) {j = 1, . . . , + 1) for 
qudit states was considered in [15| . The important property of this probability vector is that for the 
corresponding density operator p describing the quantum states the vector components occupy only a 
subset of the corresponding simplex. To see this let us consider the conditional probability 



for m = 1, . . . , and a = 1, . . . , + 1 (with N = 2j + 1). One assumes that the set of unitaries 
{Ui, . . . ,U N+i} defines a quorum, that is, one may reconstruct p out of p(m|a). Note that p{m\a) encodes 



Wim,U) = {m\UpU''\m) , 



(1) 




(2) 



p{m\a) = (m\UapUl\m) , 



(3) 
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{N + 1) X [N — 1) = N'^ — 1 real parameters, that is, exactly the same nmnber as the density operator 
p in iV-dimensional complex Hilbert space. Fixing the {N + l)-dimensional probability distribution Tr^, 
that is, we take into account the measurement of Ua with a weight tTq > 0, one defines 

pi"' :=p(fc|a)7r„ , (4) 

and finally N{N + l)-dimensional probability vector 

p=ip['\...j^\...jr'\---jr'r . (5) 

It is clear that only very special vectors P correspond to legitimate quantum states on A'^-level quantum 
system. Any such vector has to respect (N + 1) linear constraints 

pi"^ + ...+piv"^-^o , (6) 

for each a = 1, 2, . . . , + 1. It is clear that a proper set of Ua enables one to reconstruct p out of p'^"^ 
Note, that positivity of p, that is. p > 0, provides highly nontrivial constraint for a set probabilities p*-"-* 
and hence for a set of admissible vectors P. Hence, probability vectors P corresponding to legitimate 
quantum states define only a proper (iV^ — l)"dimensional convex subset Qn of (7V^ + 7V — l)~dimensional 
simplex I]jv(Ar+i). 

To summarize: using a tomographic approach with a fixed quorum Ua any quantum state p of iV-level 
quantum system may be attributed a probability vector P in N{N + l)-simplex. Vectors compatible with 
legitimate quantum states define a proper convex subspace Qn- 

As an illustration consider a qubit state p. Using the Bloch representation one has 

p = ^ (I2 + xax + yay + za^) . (7) 

The requirement p > provides the following constraint for the Bloch vector r = (a;,y, z): 

\r\^ = + <1 . (8) 



Pil\x) = p(2|x) 

pm = pim - 

p(l|z) = p{2\z): 
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where as a quorum one takes {Ux — (Jx^Uy — <Ty^Uz ~ Cz}. Hence fixing a probability distribution 
{tTx , tTj, , tTz ) one defines 

p = {p[''\pi'\p['\pi'\p['\p'^^y . (9) 

One easily finds the following relations: 

= y(1 + [Pi2 + P2i]) , P'f^ = ^{l-^[P12-P2l]) , p['^=nzpu, 



^(x) ^ _ [^^2 +P21]) , P2^^ = y(1+«[P12 -P21]) , P2^^=7r^P22 

Vectors p^,"-* satisfy 3 linear constraints 



p(-)+p(-)=^_ p(^)+p(^)=^^, (10) 
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The quadratic constraint ([S]) yields 



al al al 

where 

TTj; TTy TTz 

ax = — , ay ^ Y ' Y ■ 

Hence a Bloch ball ([5]) is transformed into an ellipsoid pT|). Note that for uniform distribution tt^, = 
TTj, = TT^ = i the above ellipsoid reduced to a ball. 



3 Classical vs. quantum dynamics of probability vectors 

Both classical and quantum states of finite levels systems can be encoded as probability vectors in the 
corresponding simplex. Now comes a crucial question: suppose we are given a probability vector Po and 
its evolution F(t) for t > 0. Can we discriminate between classical and quantum evolution? 

Let us briefly recall standard description of classical and quantum stochastic dynamics. Consider first 
an iV-levcl classical systems described by the following master equation 

^p(i) = Mpit) , p(0) = Po , (13) 

where p = (pi, . . . ,pnY denotes probability distribution, and M stands for the time-independent gener- 
ator. The solution defines a classical dynamical map T{t) such that p{t) = T{t)po. It is clear that T{t) 
is a stochastic N x N matrix such that r(0) = Ijv- Let us recall that a real n x n matrix T is stochastic 
iff: i) Tij > 0, and ii) J^i '^ij — 1 fo'" ^3 ~ I7 • ■ • : "n- The master equation rewritten in terms of 
T{t) reads 

jT[t) = MT{t), T(0) = lw, (14) 

and hence the solution is given by T{t) = e*^^. It is well known [T^ that e^^^ provides a stochastic matrix 
for alH > if and only if M satisfies the following Kolmogorov conditions : i) > for i 7^ j, and ii) 
Mij — for all J = 1, . . . , A^. In this case M possesses the following representation 

n 

M,j = TTy - (5,j TTkj , (15) 

fc=l 

with TTij > ior i ^ j. Actually, the diagonal elements ttu cancel out and hence they are completely 
arbitrary. Putting this form into (|T4)) one arrives at the classical Pauli rate equation 

d ^ 

^ft(0 = [^y^^J W " ^j^kW] > (16) 

where nij play a role of transition rates. 

Consider now the dynamics of a quantum n-level system. The corresponding quantum master equation 
reads 

^p(t) = Lp{t) , p{0) = PO , (17) 

where p denotes a density operator, and L stands for the time-independent (quantum) generator. The 
corresponding quantum evolution gives rise to the family of dynamical maps A{t) : T(H) — > T('H) such 
that p{t) = A{t)pQ. For any t > a map A{t) is completely positive and trace preserving (i.e. defines 
a quantum channel). The formal solution A{t) = e*^ defines legitimate dynamical map if and only if 



4 



the corresponding generator has a well known Gorini-Kossakowski-Sudarshan-Lindblad form [THl [IHl (see 
also HO]) 

Lp = -t[H, p] + lY. ([^fc/^' ^k] + [Vk,pV^]) , (18) 



2 

A: 

where i/^ = H <E S(H) and Vk G *B(H) are arbitrary. Interestingly, if one fixes an orthonornial basis 
{ei, . . . , e„} in H, then the following matrix 

M,, = Tr[P,L(P,-)] , (19) 
where Pi = |e,;)(ei|, satisfies Kolmogorov conditions, and 

T,,{t)^TT[P.Mt){P,)], (20) 

is stochastic for all t > 0. 

Up to now we described classical and quantum evolutions using different frameworks: stochastic 
evolution of a probability vector and completely positive evolution of a density operator. To compare 
how classical and quantum system evolve let us encode p as a n{n + l)-dimensional probability vector 
P £ Qn- It is clear that the quantum master equation ()17|) rewritten in terms of P gives rise to the 
following linear equation 

|p(t)=MP(t), (21) 

and the corresponding solution has the following form P{t) ~ T(t)Po, with T{t) — e*^. 

Now comes the question: suppose one is given a master equation (|2ip for a probability vector P G 
Qn C S]„(„-|_i) Is it possible to discriminate whether P(t) represents classical dynamics of n{n + l)-lcvel 
system or quantum dynamics of n-levcl system encoded in P? Any M satisfying Kolmogorov conditions 
gives rise to legitimate classical dynamics. It means that in the classical case T(t) ~ e*^ is always 
stochastic matrix for t > 0. Interestingly, for quantum dynamics it is no longer the case. It means 
that if T{t) represents quantum evolution it needs not be a family of stochastic maps. Equivalently, the 
corresponding generator needs not have Kolmogorov form. Hence, whenever T{t) violates stochasticity 
it proves the quantumness of the corresponding evolution. 



4 Examples: qubit dynamics 



In this section we provide simple examples of qubit dynamics giving rise to the evolution T(t) which 
needs not be stochastic. 

Proposition 1. IfLp— —i[H,p], then the corresponding M. is never of Kolmogorov form. 

Indeed, note that in this case the dynamics is invertible, that is, L and —L share the same property. 
Now, if M were of Kolmogorov form then — M is not and hence the symmetry is broken. Note, that 
M„ = tr[P,L(P,)] = 0. 

Example 1. Consider a Schrodingcr evolution of a qubit state governed by the von Neumann equation 

d 



dt 



p{t) = -i[H,p{t)] , 



(22) 



with H = Ljax ■ One finds 
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with fi ~ TTy/TTz- It is clear that M is not a Kolmogorov generator and hence the evohition T(<) = e*** is 
not stochastic. We stress that a sign of cj does not play any role. 

One might be tempted to expect that purely dissipative quantum dynamics leads to legitimate stochas- 
tic evolution T(i). To show that it is not the case let us consider an elementary generator constructed 
out of a single operator V: 

Lp = VpV'^ - ^ (y + pV^fV) . (24) 

One finds for the diagonal elements 

Pii = -7Pii +7P22 + Kpi2 + , 
P22 = IPii - JP22 - K pi2 -K P21 , 



where 



and hence 



7 = 1^12 



d (z) 
d {z) 



n = - \{'^\V^V\l) = ^(^11^;! - ^21^) 



(z) 

'7Pi 



(^) 
-IP2 



X2P[^^ 



(v) 



(z) (z) V (^) V (^) 

IPi -IP2 - Xip\ - X2P\ 



for real parameters Xi , X2 , Yi , 1^2 (they are easily calculable in terms of n and k) . It is therefore clear that 
unless Xi = X2 = Yi = I2 = the corresponding dynamics of probability vector P can not be stochastic: 
indeed if Xi > then ~Xi < and hence the generator does not satisfies Kolmogorov conditions My- > 
for i ^ j. The above analysis gives rise to the following 

Theorem 1. A legitimate quantum evolution A{t) gives rise to legitimate stochastic evolution T{t) if and 
only ifT{t) is block diagonal, that is, 

m = r(")(t) , (25) 

a 

and each T^^^t) defines stochastic evolution of the probability distribution p{m\a) ~ Pm^ /tTo- Equiva- 
lently, the corresponding generator 

M = M(") , (26) 

a 

and each M*-"-* satisfies Kolmogorov conditions. 

Interestingly, several well known quantum generators leads to legitimate Kolmogorov form of M. 
Example 2. Consider a dissipative evolution governed by the following master equation 



(27) 



where the dissipative part L^) is defined by 

LdP = ^{71 ([(7+^,(7-] + [(T-,pa+\) +72([cr-p,cr+] + [a+,pa-]) +-iz{(JzpcTz - p) } , (28) 

with cr+ ~ |2)(1| and o-_ = |1)(2| the standard raising and lowering operators. One easily find the 
corresponding generator 
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where T = 5(71 + 72) + 73 and u = iTx/TTy. Interestingly, block-diagonal blocks satisfy separately 
Kolmogorov conditions. Note that if cj = 0, i.e. dynamics is purely dissipative, then M is block-diagonal 
and hence perfectly of the Kolmogorov form. 

Example 3. Consider the following generator 

3 

Lp = ^ Jk{crkpcrk - p) , (30) 
fe=i 

with 7fc > 0. Evolution generated by the above generator is called random unitary since it is defined by 
the following formula 

3 

Ht)P = X! P^'^*^ ^^iP<^^i. : (31) 
where p^it) is a probability vector for t>{). One finds 

Po(t) = i[l+Ai2(t)+Ai3(i)+A23(0] , Pl(t) = i[l-Ai2(t)-Ai3(0+A23(t)] , 
P2{t) = i[l-Ai2(t) + Ai3(t)-A23(0] , P3(i) = ^[l + ^12(t)-Ai3(t)-A23(0] , 

where Aij{t) ~ e"^!^''''"*"''^)*! . Using the following matrix representation 

Lp = 



ill + 72)(P22 - Pll) -Pl2(7l + 72 + 273) + P2l(7l - 72) 

P2i(7i + 72 + 273) + ^12(71 - 72) (71 + 72)(pii - P22) 



one easily finds 
with 



M = Af(^) © M'^y^ © M(^) , (32) 

Af(") = ( -T- ) , A/('') ^( \ , = ( ) ' (33) 



and 



7a; = 72 + 73 , 7i; 71 + 73 , 7z = 71 + 72 



Hence random unitary dynamics of a qubit provides perfectly stochastic dynamics of the 6-dimensional 
probability vectors P. Interestingly, one has three 2-dimensional stochastic evolutions 

— p(") = M(")p("' , (34) 
dt 

with a = x,y,z, and M^") = ^aicTx ~ 12)- Note that asymptotically p^"^(t) 7rQ,(i, i)"^, i.e. each 
2-dimensional probability vector p^") /tTq, becomes maximally mixed. 



5 Conclusions 

To resume, we point out our main results. The dynamics of random classical systems is usually described 
by time dependence of a probability vector which is identified with the system state. The linear dynamics 
is associated with maps realized by stochastic matrices forming a semigroup. In the framework of the 
tomographic approach the state of quantum systems arc also identified with probability vectors. But 
these vectors are located in a subdomain of the simplex, which is the usual region for classical probability 
vectors. The quantum subdomain is determined by constraints to which obey all the probability vectors 
determining the quantum system states. The linear maps of this subdomain form a semigroup. The 
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matrices realizing the linear maps are not only stochastic matrices, as they may contain negative matrix 
elements (an analog of quantum quasi-distributions like the Wigner function). 

We analyzed quantum dynamics of the density operator in the framework of a single probability vector 
P. Interestingly, the corresponding dynamical map T{t) needs not be stochastic contrary to the classical 
evolution of P. Therefore, violation of stochasticity witnesses quantumness of evolution. It turns out 
that unitary dynamics always violates stochasticity. We showed that purely dissipative dynamics of qubit 
analyzed in Examples 2 and 3 is perfectly stochastic being a direct sum of three stochastic evolutions. It 
seems that this feature is responsible for solvability of these models. 

We will present in a future paper other examples of dynamical quantum maps belonging to semigroups 
of matrices, with both nonnegative and negative matrix elements. 
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